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ABSTRACT 

We present new 1-1.25 (z and J band) Subaru/IRCS and 2 (K ba nd) VLT/NaCo da ta for 
HR 8799 and a rereduction of the 3-5 MMT/Clio data first presented bv lHinz et al] ()2010D . Our 
VLT/NaCo data yields a detection of a fourth planet at a projected separation of ~ 15 AU - "HR 
8799e". We also report new, albeit weak detections of HR 8799b at 1.03 fim and 3.3 /im. Empirical 
comparisons to field brown dwarfs show that at least HR 8799b and HR8799c, and possibly HR 
8799d, have near-to-mid IR colors/magnitudes significantly discrepant from the L/T dwarf sequence. 
Standard cloud deck atmosphere models appropriate for brown dwarfs provide only (marginally) 
statistically meaningful fits to HR 8799b and c for unphysically small radii. Models with thicker cloud 
layers not present in brown dwarfs reproduce the planets' SEDs far more accurately and without the 
need for rescaling the planets' radii. Our preliminary modeling suggests that HR 8799b has log(g) — 
4-4.5, Te// = 900K, while HR 8799c, d, and (by inference) e have log(g) = 4-4.5, Te// = 1000-1200K. 
Combining results from planet evolution models and new dynamical stability limits implies that the 
masses of HR 8799b, c, d, and e are 6-7 Mj, 7-10 Mj, 7-10 Mj and 7-10 Mj. "Patchy" cloud 
prescriptions may provide even better fits to the data and may lower the estimated surface gravities 
and masses. Finally, contrary to some recent claims, forming the HR 8799 planets by core accretion 
is still plausible, although such systems are likely rare. 
Subject headings: 



1. INTRODUCTION 

The HR 8799 planetary sy stem is the first dire ctly im- 
aged multiplanetary system ()Marois eTan [20081) . Along 
with Fomalhaut and /3 Pic, it is also the only imaged 
system with companion mass ratios and separations rea- 
sonably close to the giant planets in the Sola r System 
(|Kalas et al.ll2008t iLaerange et al.ll2009l l2010l) ^°. After 
the initial detection of HR 8799bcd, one or more planets 
were recovered in prior datasets ( Lafrcnicrc ct al. 200^ 
Fukagawa et~al] 120091: iMetchev et al.l 1200911 . Recently, 
Marois et al.l (|2011| ) imaged a fourth planet - HR 8799e 
- which we independently detected (see §2). 

Mass estimates based on cooling models yield 5- 
11 Mj for HR 8799b and 7-13 Mj for the other 
planets ijMarois et al.l I2008L 1201 If ). Dynamical con- 
straints placed by HR 8799bcd imply that the com- 
panion s likely have masses below the deuterium-burning 
limit (jSpiegel et al.l l2010D and are kept stable by 
resonant interactions (iFabrvckv and Murrav-Clavl 120101 : 
iMoro-Martin et al.l |2010D . Including the fourth planet. 
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iMarois et all (|201lD argue that the planets most likely 
have masses at the low end of the range allowed by cool- 
ing models. With masses of « 5-13 Mj, the HR 8799 
planets then bridge the gap between the solar system's 
gas giants/ Jupiter- mass pla nets detecte d by radial ve- 
locity surveys (e.g. Howar d et al.l l2010l ) and low-mass 
brown dwarf companions to nearby stars such as GJ 758B 
and PZ Tel (iThalmann et al.l l2009l : ICurrie et all 120101: 
iBiller et al.ll20T?ii r 

Recent studies complicate our understanding of the re- 
lationship between brown dwarfs, the gas giants detected 
in RV surveys, and the HR 8799 planets. The planets' 
masses are significantly larger than most pla nets detected 
by ra dial velocity and transit methods. iMarois et al.l 
(|2008() noted that the planets appear slightly redder than 
the distribution of H/H-K^ colors for old field brown 
dwarfs. The K-band spectrum of HR 8799b is not well 
matched by typical L and T-type brown dwarf spectra 
(jBowler et al.ll2010D . 

Comparisons between the HR 8799 planet photom- 
etry /spectroscopy and atmosphere models reveal ad- 
ditional difficulties in understanding their properties 
within the theoretical framework of standard, cloud deck 
models that trac k the field L/T dwa. rf sequence. In the 
discovery paper, IMarois et al.l (|2008[ ) briefiy mention a 
discrepancy between temperatures derived from atmo- 
sphere models and those estimated from more simple, 
and presumabl y most accurat e , cool ing model estimates. 
More recently, I Bowler et al.l (|20100 provide a detailed 
comparison between the HR 8799b spectra and 1.1-4.1 
/im photometry and predictions from standard atmo- 
sphere models. They show that the 'best-fit' tempera- 
tures derived from modeling are inconsistent with cool- 
ing model estimates. They also explicitly show that the 
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implied radii for best-fit models are well below the 1.1- 
1.3 Rj range allowed by standard cooling models (e.g. 
0.3-0.5 Rj). 

To interpret these modeling difficulties, iBowler et al.l 
()2010l ) argue that a different atmospheric structure, 
namely atmospheres with stronger cloud coverage, may 
better explain the HR 8799b SED. Since atmospheric 
dust entrained in clouds absorbs more efficiently at 
shorter wavelengths, photometry for HR 8799b at wave- 
lengths shortward of J band would provid e a crucial test 
of th e planet 's level of cloud cove rage (cf. iBurrows et al.l 
I2006f l. The IBowler et all ()2010f ) study also found dif- 
ficult y in reconcil i ng th eir model fitting of detections 



fromlMarois et al.l (|2008[ ) with 3-5 fj,m upper limits from 



iHinz et al.l ( 201(3 ). More sensitive photometry at these 



wavelengths would then provide better modeling con- 
straints. 

In this study, we investigate the atmospheres and dy- 
namics of the HR 8799 planets using new observations 
obtained at the Subaru Telescope an d VLT and a rere - 
duction of MMT data presented by iHinz et al l (|2010D 
Comb ined with photometry presented bv iMarois et al.l 
(|20f)8h . our data yield nine photometric points spanning 
1-5 fim for a detailed comparison to the properties of 
field brown dwarfs. This wavelength range also provides 
a sensitive probe of the effects of surface gravity, temper- 
ature, (non) equilibrium chemistry, metallicity, and cloud 
coverage. 

We compare the planets' SEDs to atmosphere mod- 
els exploring a phase space defined by these effects. By 
quantifying the model fits, we determine the range of pa- 
rameter space that fails to characterize the planets' SEDs 
and identify the subset of models that more accurately 
reproduce the data and may better represent their atmo- 
spheres' physical properties. These results will then be 
used to more thorougly and accurately probe the plan- 
ets' atmospheric properties in a companion paper (Mad- 
husudhan et al. 2011, in prep.). 

Our study is structured as follows. §2 describes our 
observations, image processing, and detections for each 
dataset. The first part of §3 compares the HR 8799 
planet photometry to the L/T dwarf sequence and the 
IR properties of other very low-mass objects (M < 25 
Mj). The rest of §3 presents preliminary comparisons 
between the HR 8799 planet SEDs and planetary atmo- 
sphere models. §4 describes simple dynamical modeling 
of the system to identify the range of masses for dynam- 
ically stable orbits. §5 summarizes our results, discusses 
our work within the context of previous studies of HR 
8799 and planet imaging in general, discusses how our 
results fit within the context of planet evolution models, 
and comments on the plausible formation mechanism(s) 
for the planets. 

2. DATA 

2.1. Observations 

Our study combines data from three facilities - 
VLT/NaCo, Subaru/IRCS, and MMT/Clio - at six 
broadband filters centered on 1.03 fim to 4.8 iim. The 
VLT data are the most sensitive and were obtained to 
place limits on the presence of other candidate planets 
in the system. The Subaru data at 1.03-1.25 fim were 
taken to probe the effect of clouds on the planets' atmo- 
spheres. Finally, we rereduced the 3-5 fj,m MMT/Clio 



data first presented bv lHinz "eraII((20Tl using our reduc- 
tion pipline, which utilizes advanced image registration, 
PSF removal, speck le suppression, and flrp c calibration 
routines (e.g. LOCI iLafreniere et al.ll2007al) also used in 
IMarois et all (jl008). 

All of our data were taken in angular differential imag- 
ing mode (jMarois et al.l l2006l) . where the instrument 
rotator is adjusted to stay at a fixed angle with re- 
spect to the (changing) parallactic angle, resulting in 
the field of vi ew rotating with time. Combined with the 
IMarois et al.l data, we thus have data spanning nine pho- 
tometric filters that is largely reduced self consistently. 
Table [1] summarize basic properties of our observations. 

2.1.1. VLT/NaCo K, hand Data 

HR 8799 was imaged with VLT/NaCo on six separate 
nights in October 2009 as a part of a separate study of the 
HR 8799 planets (RI. Daniel Apai; Apai et al. 2011, in 
prep.). Once publically available, the science and calibra- 
tion data were downloaded from the ESO VLT archive 
for October 8-11, nights over which the field rotation for 
the HR 8799 data was > 30-45 degrees, resulting in a 
small (r ~ 0.22") inner working angle. 

The data were taken with the 13.27 mas pixel scale, 
without coronographic masks, and in pupil tracking 
mode allowing for angular differential imaging. All data 
consist of coadded 0.345s exposures totaling ~ 43s a piece 
and are stored in the standard NaCo datacube format. In 
this paper, we focus specifically on the October 8 data, 
which had the highest quality and greatest amount of 
field rotation. Apai et al. (2011, in prep.) will later 
present a larger study combining all October 2009 data 
and FaU 2010 data. 

2.1.2. Subaru/IRCS z(Y) and J band data 

HR 8799 was targeted for direct imaging on August 
15, 2009 with the Subaru Telesc ope using the Infrared 
Camera and Spectrograph flRCS^ lTokunaga et al.|[r998D 
and AO-188 adaptive optics system in natural guide 
star mode. The data were taken in the Mauna Kea 
J band filter (~ 1.25 /xm) and the z filter centered on 
1.033 /i m, analogous to t he better- known Y band fil- 
ter fe.g. iHillenbrand et al.l l2b02l" During our observa- 
tions, conditions were photometric with fair natural see- 
ing 0.65-0.75"). AO-188 yielded a corrected image 
with FWHM(PSF) ~ 0.06" in z and 0.064" in J. For ah 
observations, the native pixel scale is 20.57 mas/pixel; we 
used the 0.8" diameter, non-transmissive coronographic 
mask to block most of the primary starlight. 

The z data were taken using 30 second exposures con- 
sisting of 6 coadded frames to avoid saturation at separa- 
tions corresponding to HR8799bcd for a total integration 
time of 4500s. The J band data consist of 25 second coad- 
ded exposures for a total integration time of 1080s. The 
z data were observed through transit yielding a total field 
rotation of 172 degrees. The J data were taken about an 

The zeropoint wavelength for the z filter listed on the IRCS 
webpage is 1.033 /im with a width of 0.073 ^tm. The Y-band filters 
for comparable cameras are slightly wider but otherwise quite sim- 
ilar: filters for Keck/NIRC (there called "Z"), UKIRT/WFCAM 
and Gemini/NIRI have zeropoint wavelengths of 1.032, 1.031 and 
1.02 /im and widths of 0.156, 0.1 and 0.1 /im, respectively. The 
IRCS z filter should not be confused with the Sloan z' filter, which 
covers shorter wavelengths. 
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hour after transit resulting in very poor field rotation (~ 
6.4 degrees). 

2.1.3. MMT/Clio 3-5 fim Data 

MMT/C l io obs ervations were previously described by 
iHinz et all (|2010^. Briefly, HR 8799 was imaged in three 
separate observing runs - November 21, 2008, January 
10, 2009, and September 12, 2009 - at the L' (3.8 fj.m) 
and Barr M (4.8 fim) and a shorter wavelength filter cen- 
tered on 3.3 iim methane absorption feature and extend- 
ing from 3.12 fim to 3.53 fim. We focus on the November 
2008 and September 2009 runs, which had sufficient field 
rotation for angular differential imaging. The pixel scale 
for all Clio data is 48.57 mas/pixel. 

The [3.3], L', and M data were imaged for 6780s, 5690s, 
and 9600s: the total field rotation for data in these three 
filters was 125.3 degrees, 72 degrees, and 31.8 degrees. 
While observing conditions for the L' and M data were 
clear, the [3.3] micron data were taken through variable 
seeing in two sets between which the AO system failed 
to yield an acceptable correction. 

2.2. Image Processing/Data Reduction 
2.2.1. Basic Image Processing and Image Registration 

For our Subaru/IRCS and VLT/NaCo data, we first 
performed standard dark subtraction, fiat fielding, and 
bad pixel masking. While the NaCo data followed a four- 
point dither pattern which should wash out image dis- 
tortion errors, the IRCS data were not dithered. Thus, 
each IRCS frame was corrected for image distortion us- 
ing polynomial fits, resulting in a revised pixel scale of 
20.53 mas/pixel. 

For the MMT/Clio data we then performed sky sub- 
traction. We constructed Clio sky frames from median- 
combined images obtained for each nod position and sub- 
tracted to remove the sky background. Final pixel values 
for each VLT/NaCo image were nominally constructed 
from the average pixel value drawn from each frame in 
the datacube. For regions within 1" of estimated stellar 
centroid, we determined the average pixel value after it- 
eratively clipping 5a outliers. For all datasets, bad pixels 
were identified as outliers within a moving-box median 
filter, fiagged, and interpolated over. 

Our image registratio n proce dure closely f ollows that 
of iLafreniere et all (|2007b^ and .Marois et al.l (j2008,). We 
first copied each image into a larger blank one, coarsely 
registering them using a priori knowledge about the cen- 
ter of the coronographic mask (for IRCS) or a gaussian fit 
to a convolved version of the image using the IDL func- 
tion gcntrd.pro. For precise image registration, we center 
one image using a 2D cross correlation function relating 
it to a 180 degree rotation of itself. We then identify the 
fractional pixel offsets between the reference image and 
subsequent images yielding the highest correlation. The 
region of interest used to register IRCS images is focused 
on diffracted light from the secondary spider. For the 
Clio and NaCo images we used the non-saturated por- 
tions of the stellar PSF, since the diffracted light from 
the spider is highly suppressed. 

2.2.2. Localized Combination of Images (LOCI) Speckle 
Supression Processing 



Further data reduction f ollows the ADI/LOCI re duc- 
tion procedur e descr ibed bv lLafreniere et al\ (|2007al ) and 
iMarois et al.l ()2008[ ). We first subtract out the time- 
independent component of the stellar PSF, exploring two 
methods. In the first method, we median combine all im- 
ages for a reference PSF which we sub tract from each im - 
age, the simple ADI method used by iHinz et all ()2010D . 
In the second method, we construct a two-dimensional 
radial profile for each image and subtract it to remove 
the smooth seeing halo. 

Next, we perform the LOCI speckle suppression algo- 
rithm on the residual images, derotate the processed im- 
ages and median combine them for a final science image. 
We compared reductions for a range of LOCI input pa- 
ramet ers - dr, N5, Na, and geom (see ILafreniere et al.l 
l2007aL for definitions) - to identify the set that max- 
imized the signal -to-noise of the planets, using the set 
recommended by ILafreniere et al.l (|2007al) as a starting 
point. Our pipeline also produces the simple ADI reduc- 
tion as a byproduct, useful for a separate, sensitive iden- 
tification of HR 8799b, whose detection in some filters 
(e.g. [3.3], M) may be more severely limited by photon 
noise than by speckle noise. 

2.3. Planet Detections and Astrometry 

To identify detected planets in our images, we com- 
pute the standard deviation and si gnal-to- noise ratio 
of pixel values in con centric annuli ( Currie et al.ll20iol : 
iThalmann et al.ll2009l ). As a check on our results, we 
compare our astrometry for candidate detections in a 
give n filter with that obtained by us in other filters and 
from IMarois et al.l (|2008l 1201 iD during the Fall 2008 and 
2009 epochs. We claim a detection of a planet indepen- 
dent of other datasets ff SNR > 5. For SNR = 3-5, the 
centroid of the candidate planet detection must be the 
same as that reported for the planet data where SNR 
> 5 within astrometric errors (typically 0.5 pixels). We 
centroid the planet using the IDL functions gcntrd.pro 
and cntrd.pro and adopt a minimum astrometric uncer- 
tainty of 0.5 pixels to account for image registration and 
centroiding/orientation errors. The rightmost column of 
Table [T] summarizes our planet detections and Table [2] 
lists their astrometry. 

2.3.1. VLT/NaCo Detections 

Figure [1] shows our reduced VLT/NaCo Kg band im- 
age. HR 8799 b and c are detected at better than 25 
a, while HR 8799d is detected at 10 a. The planets are 
also free of deep, negative flux troughs at the same sep- 
aration but slightly different position angles that results 
from LOCI being applied to datasets with poor field rota- 
tion or those where most exposures are take n well before 
or well after transit (e.g. IMarois et al.ll20l"ol) . 

Additionally, our data show a detection of an another 
point source located interior to and in the same quad- 
rant as HR 8799d consistent wi t h a fo urth planet - "HR 
8799e" . Recently, IMarois et al] ()201lD announced a mul- 
tiepoch detection of HR 8799e in K-band and L'-band 
using Keck/NIRC2. Their detection significance in K- 
band using Keck/NIRC2 is slightly better than ours (^ 
5 (7 vs. our ~ 4 cr). Our photometry using methods de- 
scribed in §2.4 yields an absol ute magnitud e of m{Ks) — 
12.89 ± 0.26, consistent with IMarois et all 's estimate of 
12.93 ± 0.22. 
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Figure [2] compares our astrometry. We measure a cen- 
troid position of [E,N] = [-0.306" ± 0.007", -0.217" ± 
0.007"], implying a projected separation of 14.8 AU ± 0.4 
AU. The a verage of the August and November 2009 posi- 
tions from lMarois et al.l ()201lD is [-0.304,-0.203] with an 
intrinsic uncertainty ~ 0.01" . Our position is then con- 
sistent with theirs to within 1.4 a. Ou r imphed proje cted 
physical separation is consistent with iMarois et al.l 's es- 
timates from multiepoch data: 14.5 AU ± 0.5 AU. 

2.3.2. Suharu/IRCS Detections 

Figure [3] show reduced images at J and z obtained with 
IRCS. In spite of poor field rotation severely limiting the 
performance of the ADI/LOCI processing and precluding 
detectability of objects within ^ 1", we clearly detect 
the b planet in our J band data at a ~ 10 ti significance 
(top panels). In spite of good seeing conditions, good 
field rotation, and a 70-minutc integration time, we fail 
to detect any of the planets at a > 5 cr significance in z 
band (bottom panels) . Our reduced image reveals a weak 
detection of HR 8799b with SNR - 3.7 (bottom-right 
panel) and a centroid within 0.25 pixels of its centroid in 
the J-band data obtained one hour later with the same 
instrument. However, we do not detect HR 8799c or 
d in our z data. To verify that our low signal-to-noise 
detection of HR 8799b and nondetections for the other 
planets do not result from errors in derotation or a jump 
in parallactic angle near transit^^, we introduced fake 
planets into each image with a flux equal to ~ 10 times 
the local noise of the final image, reran our reduction 
pipeline separately for frames before and after transit, 
and recovered their detections. 

2.3.3. MMT/Clio Detections 

Figures m [51 and H] show reduced images obtained with 
MMT/Clio in the L', [3.3], and M fihers. In the L' fil- 
ter, we detect HR 8 799bcd with signa l-to- noise higher 
than that reported bv lHinz et al.l (|2010[ ) and comparable 
to that obta i ned in shorter Keck/NIRC2 exposures by 
IMarois et all (|2008( ). In the [3.3] filter, we detect the c 
planet at SNR > 5. We m arginally detect the b planet 
at SNR 3.8. iHinz et al.l (|2010) formally report a non- 
detection for b at [3.3] as they adopt a 3cr threshold for 
detections, though they identify a cluster of pixels 2.8 
a above the background consistent with b and roughly 
coincident with our MMT/Clio and VLT/NaCo detec- 
tions. 

Converse l y, we do not detect HR 8799d in [3.3], while 
iHinz et all (|2010[ ) report a low-significance detection. 
This disagreement is surprising since LOCI greatly im- 
proves the planet sensitivity at small separations such as 
that for the d planet (Lafrcnicrc et al. 2007a). Further- 
more, there is a 40mas offset be tween the reported HR 
8799d centroid from IHinz et al.l and that from our 10 a 
VLT/NaCo detection obtained three weeks later. While 
their detection is likely instead residual speckle noise, our 
qualitative conclusion that HR 8799d is very fain t at 3. 3 
/im is consistent with theirs. As with lHin"z eTall (l2010h . 
we do not detect any of the planets at M band. 

For an example of this phenomenon, see 
|http: //www2.keck.hawaii.edu/inst/nirc2/vertAngJump.html| 



2.4. Photometry for Detections and Upper Limits for 
N on- Detections 

Photometry for each dataset was performed with 
IDLPHOT with the aperture radius set to the 
0.5xFWHMi,„age- In aU exposures, the stellar PSF 
core is either saturated or obscured. For initial pho- 
tometric calibration, we obtained observations of the 
stellar primary viewed through a neutral density filter 
(MMT/Clio, VLT/NaCo) or observed standard stars im- 
mediately prior to and after our science exposures (Sub- 
aru/IRCS). 

Faint companions to stars observed in ADI and pro- 
cessed with LOCI lose flux due to field rotation and self- 
subtraction. To further calibrate our photometry, we 
introduce and measure the flux for faint point sources at 
random position angles over separations encompassing to 
the HR 8799 planets (0.25"-2") in each registered frame, 
rerun our ADI/LOCI pipeline, compute the attenuated 
flux in the final, processed images, and correct for this 
attenuation. Figure [7] illustrates this flux loss, compar- 
ing the input and output flux for fake points sources for 
our MMT/Clio L' data. While images processed using 
a simple ADI reduction lose ~ 20% of their flux, self 
subtraction is stronger with LOCI, especially at separa- 
tions less than 0.75". The attenuation curves obtained 
for data in other fllters do not differ qualitatively: LOCI 
always attenuates more flux than a simple ADI reduc- 
tion and attenuation is signiflcantly more severe at small 
separations. 

To place limits on our nondetections, we compute 
3(7 upper limits where we correct our nominal sensitiv- 
ity limits for point source self-subtraction inherent in 
ADI/LOCI. The noise is defined in concentric annuli as 
before, since in most cases (for HR 8799 c and d) radially- 
dependent speckle noise dominates over radially indepen- 
dent photon noise. Despite using LOCI, our detection 
upper limit at 3.3 jim f or HR 879 9 d is b righter than 
the magnitude listed by iHinz et al.l ()2010[ ). Moreover, 
our upper limits for HR 8799b cd at M are consistently 
brighter than those reported by iHinz et ahl . 

In both cases, the disagreement is likely explainable 
by our correction for point source self subtraction in de- 
riving upper limits from the standard deviation of pixel 
values. iHinz et al.l construct a reference PSF by median- 
combining all frames and then subtract this reference 
PSF from each image. For the 3.3 /im data, our reduc- 
tion pipeline predicts that this processing should atten- 
uate about half of the point source flux at HR 8799d's 
position^-^ . For the M band data, fleld rotation is poorer 
and thus self subtraction with this reduction procedure 
is severe, reaching over 75% at the position of HR 8799d 
as nearly half the frames are obtained ~ 3 hours after 
transit and thus at essentially one position angle. Thus, 
the gain in sensitivity due to LOCI is reduced by self 
subtraction, resulting in brighter 3 cr upper limits. 

3. PHOTOMETRIC ANALYSIS: CONSTRAINING THE 
ATMOSPHERIC PROPERTIES OF THE HR 8799 
PLANETS 

Com bining our data with that from IMarois et al.l 
(|2008l) yields planet flux measurements at nine separate 

^•^ While the total field rotation is large, ~ 127 degrees, the vast 
majority of the frames were taken over a time interval mth only 
30 degrees of field rotation. 
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wavelengths from 1 to 5 fj,m. In this section, we use this 
rich muhiwavelength samphng of HR 8799's planet SEDs 
to provide an empirical comparison with other cool, 
substellar-mass objects and simple atmospheric model- 
ing constraints on the planets' properties. 

3.1. Near-to-Mid IR Colors of the HR 8799 Planets 
3.1.1. Methodology 

To compare the near-to-mid IR properties of the HR 
8799 planets with those for other cool, substellar objects, 
we primarily use the sample of L/T dwarfs compiled by 
iLeggett et all (|2010D The L/T dwarf sequence defined 
by the ILeggett et al.l sample allows us to determine how 
the HR 8799 planet SEDs deviate from those for brown 
dwarfs of similar temperatures. To explore how the HR 
8799 planet SEDs compare to those for other planet-mass 
objects and very low-mass brown dwarfs with Te/ / — 
800-1800K, we include 2M 1207b (5 Mj), 1RXJ1609.1- 
210524 (9 Mj), AB Pic (13.5 Mj), and HD 2030301 ) 
25 M.j) (iChauvin et aLl[200l iLafreniere et all [2?)l?)l : 
IChauvin e t al.l 120051 : iMetchev and Hillenbrandl l2006f ). 
Table m lists photometry for these objects. 

We use color-magnitude diagrams constructed from the 
Y, J, H, Ks, and L' filters to determine whether the HR 
8799 planets are similar to or under/overluminous com- 
pared to the Leggett L/T dwarf sequence. For simplicity 
and because there is no published response function for 
the IRCS z band filter, we treat the IRCS z-band mag- 
nitudes/upper limits for the HR 8799 planets as synony- 
mous with its Y-band magnitude. 

To provide a physical point of reference for the L/T 
dwarf sequence and the HR 8799 color-magnitude po- 
sitions, we overplot loci for sta ndard, chemical equ i- 
librium atmosphere models from iBurrows et al.l (|20Q6l ). 
We specifically choose the Model E case, which assumes 
that the clouds are confined to a thin layer, where the 
thickness of the fiat part of the cloud encompasses the 
condensation points of different species with different 
temperature-pressure intercept points. Above and below 
the flat portion, the cloud shape function decays to the -6 
and -10 power. Thus, above and below the flat portion, 
the clouds h ave scale height s ^ l/ 7th and 1/1 1th that of 
the gas. See IBurrows et all (|2006D for more details. 

3.1.2. Results 

Figure [8] shows our color comparisons. At least three of 
the HR 8799 planets have Kg/Kg-L' positions (upper-left 
panel) roughly consistent with those for the Leggett L/T 
dwarf sequence and with 2M 1207b. HR 8799cde have 
positions overlapping with objects near the L/T dwarf 
boundary; HR 8799b has a similar K^-L' color but is 
underluminous compared to the three other companions 
and 2M 1207b by a factor of two. It is unclear how its 
position compares to those for fleld L/T dwarfs because 
the sequence is poorly sampled at HR 8799b's band 
magnitude. 

The other three panels of Figure [8] clearly show that 
HR 8799c, d, and especially b have near-IR colors that 
depart from the L/T dwarf sequence. In J/J-H and J/J- 
K^, the L/T dwarf sequence turns towards blue colors by 
up to 1.5 mag starting at the L/T transition. While HR 
8799c's position is roughly coincident with TO dwarfs, 
HR 8799b and d follow an extension of the slope of the 



L dwarf sequence between J/[J-H,K3] = 11/[0.6,1.2] and 
15/[1.2,2] towards fainter magnitudes and redder colors. 
HR 8799d's position coincides with that of HD 203030b, 
while HR 8799b is located closest to 2M 1207b. The 
H/Y-H color- magnitude diagram shows that HR 8799c 
also is likely red/ underluminous; HR 8799b is 2.5 magni- 
tudes too red in Y-H for its H-band magnitude, indicat- 
ing that it is underluminous compared to the L/T dwarf 
sequence at both Y and J. 

Figure [5] overplots loci of standard lBurrows et al.l mod- 
els for parameters covering a range expected for low- 
mass, cool brown dwarfs - Te// = 800-1800K, log(g) = 
4-5 - and two metallicities (solar and 3x solar) With 
the exception of some L/T dwarf transition objects, the 
dispersion in color-magnitude positions for L/T dwarfs 
is well reproduced by model atmosphere loci. This in- 
dicates that L/T dwarf atmospheres can be explained 
within the phase space encompassed by the models' as- 
sumed cloud str ucture and ran^e in te mperature, gravity, 
and mctallicity (jBurrows et a l. 2006) . 

The HR 8799 planets, especially HR 8799b, are differ- 
ent. They consistently lie below the region enclosed by 
the standard model atmosphere loci, indicating that their 
near-IR luminosities are weaker compared to luminosities 
expected if their cloud structure were well represented by 
the models. HR 8799b in particular probes a completely 
different range of parameter space, lying 0.75 magnitudes 
or more redder than any standard atmosphere predic- 
tion regardless of temperature. Thus, Figure [8] suggests 
a strong contrast between the atmospheric properties of 
L/T dwarfs and the HR 8799 planets. 

To summarize, all three HR 8799 planets - especially 
HR 8799 b - have near-IR colors that cannot be easily 
understood within the field L/T dwarf sequence. They 
are consistently red and underluminous at Y and J, indi- 
cating that the 1-1.25 /im portion of their SEDs are sup- 
pressed in fiux. The HR 8799 planets also lie well outside 
the loci of standard atmosphere models used to interpret 
the physical properties of L/T dwarfs. Thus, the HR 
8799 planet atmospheres are not simply 'scaled down' (in 
mass) versions of the atmospheres of field brown dwarfs 
defining the L/T dwarf sequence. 

On the other hand, the planets' atmospheres show 
strong similarity to those for planetary-mass/low-mass 
brown dwarf companions to nearby stars. Specifically, 
HR 8799 c and d have similar near-IR colors to HD 
203030b, while HR 8799b consistently shows near-IR col- 
ors similar to 2M 1207b. The planetary-mass compan- 
ions lRXJ1609.1-210524b and AB Pic b are also red- 
der in near-IR colors compared to the L/T dwarf se- 
quence but not underluminous. Within the narrow con- 
text of our analysis, planetary-mass companions in gen- 
eral might not follow the L/T dwarf sequence. 

3.2. Fiducial Model Atmosphere Fits to the HR 8799 
Planet SEDs 

Our color comparisons motivate a further investiga- 
tion of the HR 8799 planet SEDs to better understand 
the source of the differences between their near-IR col- 
ors and those for field L/T dwarfs. To further explore 

^* We include the 3x solar models because they produce redder 
near-IR colors and the HR 8799 planets are red in the near-IR 
compared to the L/T dwarf sequence. 
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the physical properties of the HR 8799 planets we com- 
pare their photometry to atmospheric models. Because 
the color-magnitude comparisons indicate that standard 
model atmospheres provide poor fits to the planet data, 
we introduce a new set of models to explore additional 
phase space not covered by the standard models, specif- 
ically a different cloud structure: 

• The iBurrows et al.l (^200 &) Model A Thick 
Cloud Layer prescription - Like the Model 
E case, this model defines a cloud base at the 
high temperature interception point with the shape 
function at higher temperatures/pressures decay- 
ing to the -10 power. However, the cloud den- 
sity tracks the gas density at lower tempera- 
tures/pressures (si = in their terminology). 
Thus, clouds in this case are far more extended high 
in the atmosphere than in the standard Model E 
case. 

As noted in lBurrows et al.l (|2006f ). these models are 
qua litatively similar to the AMES-DUSTY mod- 
els (jAllard et al.ll200ll ). However, they are bluer 
and bri ghter than AMES-D USTY in the near IR 
because fAUard et"al] (|2001h adopts the interstellar 
medium particle size distribution. The Model A 
case fails to reproduce the L/T dwarf sequence as 
it is consistently too red a nd underluminous in IR 
color-magnitude diagrams (jBurrows et al.ll2006[ ). If 
the HR 8799 planets have thicker clouds than L/T 
type brown dwarfs, these models - or some hybrid 
between them and the " E" models - should repro- 
duce the planets' SEDs far better than the Model 
"E" case alone. 

Changing the cloud prescription radically alters the en- 
tire shape of the SED (Figure 0). The K and L' band 
fluxes are similar. However, the Model A/thick cloud 
prescription is underluminous over the Y and J pass- 
bands by an order of magnitude, underluminous at 1.65 
/zm by a factor of two but overluminous in the 3.3 fj,m 
region covering the trough produced by methane absorp- 
tion in the Model E cloud prescription. Overall, the 
Model A SED is much flatter from 1 to 4 /im. Addition- 
ally, the Model A prescription washes out the methane 
absorption feature at 1.65 /zm used to id entify the L/T 
dwarf transition (see also discussion in IBurrows et al.l 
12006) . 

Both the standard models an d thick cloud la.yer mo dels 
use the formalism described in lBurrows et al.r(|2006 ^) for 
temperatures Teff = 700-1800 K, gravities with log(g) 
= 3.75-5, and solar/super-solar abundances of metals. 
For both models, we assume modal particle sizes of 60 
/zm-100 u rn and a particle siz e distribution appropriate 
for clouds (lDeirmendiianlll964D. For bot h models we also 
assume radii from IBurrows et all (|1997D . 

3.2.1. Fitting Method 

Our atmosphere model fitting follows a s implified ver- 
sion of the fitting procedure employed by iBowler et al.l 
dJolO) to model the near-IR spectrum and photometry 
for HR 8799b. Nominally, we quantify the model fits 
with the statistic, 

n 

X — y ] {fdata,i ~ Fmodel.i) /'^data,i- (1) 
i=0 



We weight each datapoint equally. To account for vari- 
ability in emission and absolute calibration uncertain- 
ties, we set a 0.1 mag floor to a for each datapoint (see 
iRobitaille et al.ll2007l) . Because of incomplete line lists 
near the 1.6 /im CH4 band, we do no t compare the mod- 
els to data at the CH4I filter (see iBowler et al.l 120101 : 
ISaumon et"all 120071: iLeggett et all I2007D . However, we 
confirmed that this choice has no consequential bearing 
on our results. 

The z, [3.3], and M photometry include many non- 
detections. We quantitatively incorporate nondetections 
in the following way. For model predictions consistent 
with the 3(7 upper limits estimated for each nondetec- 
tion, we treat the model as perfectly consistent with the 
data and do not penalize the value. For model pre- 
dictions inconsistent with the 3ct upper limits, we do not 
automatically discard the model. Rather, we penalize 
the value by determining the flux ratio between the 
model prediction and the 3 a upper limit. Specificially, 
a model prediction 2 and 4 times brighter than the 3 a 
upper limit will be contribute 12 (3x4) and 48 (3x16) 
to the final value, respectively. 

We fit atmosphere models in two cases. First, to pro- 
vide a straightforward comparison between our data and 
the luminosity and colors predicted from atmosphere 
models we keep the radii fixed to the Burrows et al] 
(|1997( ) dwarf radii. Second, we vary the radius and iden- 
tify the scaling factor, Ck = (Rscaied / Rnominai) , that 
minimizes for a particular model: 

n 

f data,iFmodel^i / ^ florid j 

Cl = '-. (2) 

Ep2 12 
modelAI dataA 

i=0 

We nominally on ly allow the radius t o vary by ± 10% 
from the assumed IBurrows et al.l ()1997l ) values to encom- 
pass the range of radii for 5-20 M j objects at 30-300 Myr 
(~ 1.1-1.3 Rj). 

We determine which models are formally consistent 
with the data by comparing the resulting x^ value to 
that identifying the 3 and 5 a confidence limits. For the 
first case, where the planet radius is fixed, the appropri- 
ate x^ limits are 21.85 and 41.80 for 8 datapoints and 
seven degrees of freedom. For the second case - a vari- 
able planet radius - the limits are 20.1 and 39.4 for 8 
datapoints and 6 degrees of freedom. 

To select the best-fit models, we follow IBowler et al.l 
(j201Cl() by identifying the model with the smallest x^ 
and computing the Ax^ limit for a 3 ct confidence limit. 
'Best-fit' models satisfy Xmoderxtest ^99. 73%- We do 
this separately for the Model A and E cloud prescrip- 
tions. 

3.2.2. Results for Standard Cloud-Deck Models 

Table |S| summarizes our entire fitting results for models 
with the standard cloud deck prescription. Figure [10] 
displays som e of these fitting results with the planet radii 
fixed to the IBurrows et al.l (|1997l ) values. The top-left 
panel shows the distribution of x values for HR 8799b; 
the top-right panel compares the HR 8799b SED to the 
'best-fit' model. 

For each planet, the models with the lowest x^ values 
have temperatures within lOOK of those derived from 
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cooling models: T^f f = 900K, 12qOK, a nd HOOK for HR 
8799b, c, and d (see lMarois et alJl2008l ). Models with a 
3x solar abundance of metals have marginally smaller 
values. Adopting the Ax^ criterion, xll^n + Xgg.rso/^, 
the minimum values for modeling b, c, and d are 
300.9, 133.2, and 38.9. The range of temperatures and 
gravities fulfilling this criterion are Te// = 900-lOOOK, 
1100-1300K, 1000-1300K and log(g) = 4.5-5, 4.5-5, and 
4-5 for the b, c, and d planets. 

However, the fits are quantitatively very poor for HR 
8799c and (especially) b. As shown by Figure[TU] (top- left 
panel), the minimum x^ value for HR 8799b is a factor 
of ~ 5.5 times higher than the formal 5 cr confidence 
limit. The minimum x^ value for HR 8799c is twice as 
large. The large x^ difference between that for 'best-fit' 
models and the formal 5-cr confidence limit suggests that 
the models do not provide meaningful fits to the data. 
Fits to the HR 8799d SED are not quite as poor but 
include only one model with x^ < Xgg 73% • Allowing the 
planetary radii to vary over the range plausible for 5-20 
Mj objects does not qualitatively improve the model fits 
for the b and c planets (Table [5]) . 

The top righthand panel and lower panels of Figure 
[TOl illustrate how the models fail to reproduce the SEDs 
of HR 8799bcd. For example, for HR 8799b the 'best 
fit' model provides a good estimate of its Kg band and 
L' band fluxes and is consistent with its upper limit at 
M band. At [3.3], however, the model predicts too deep 
of a trough due to methane absorption, underpredicting 
the flux by a factor of ~ 3-4. Most strikingly, the model 
overpredicts the flux at Y and J band by over an order 
of magnitude. The model overestimates the H band and 
CH4S flux by a factor of ~ 2. Compared to the best- 
fitting models, HR 8799c also has too strong of a 3.3 fim 
flux and too low of a Y band upper limit. 

For modeling results discussed in Figure [TOl the scal- 
ing factors for the radii are almost always Ck =0.9 for 
temperatures greater than those predicted from cooling 
models and 1.1 for lower temperatures. To see which 
radii formally yield the smallest value s, we allow the 
radius to vary between 0.2 and 2 times the lBurrows et all 
()1997() values. The resulting trend of vs. Te// for all 
planets changes, as the minima are systematically pushed 
towards higher T^// (Te// = 1300-1400 K). However, ra- 
dius scale factors for the best-fit models imply that the 
planets are unphysically small ~Rb c d ^ 0.4, 0.6, and 0.7 
Rj. 

In summary, atmosphere models with standard, cloud 
deck prescriptions appropriate for brown dwarfs only pro- 
vide statistically meaningful fits to HR 8799b and c for 
unrealistically small radii (see also Bo wler et al.l 120101 
for HR 8799b). Assuming radii characteristic of planet- 
mass objects, we fail to find a single model that pro- 
vides a statistically meaningful fit to the HR 8799b and c 
data indicating that such models provide a poor descrip- 
tion of the planets' atmospheres (see als o iMarois et al.l 
l2008HJanson et al.ll2010HHinz et al.ll2010[ ). These resuhs 
are independent of surface gravity for log(g) = 4-5 and 
whether the planets have solar or 3x solar metallicity. 
These results then motivate us to see if models with dif- 
ferent cloud prescriptions fare better in reproducing the 
SEDs of HR 8799bcd. 



3.2.3. Results for Thick Cloud Layer Models 

Figure [TT] shows and Table [5] summarizes our fitting 
results for the thick cloud layer models. Best-fit models 
for the HR 8799 planets cover a similar range in Te// 
as the standard model fits and cooling model predic- 
tions. For HR 8799b, the best-fit model assumes Te// 
— 900K and log(g)=4.25; the range of best-fit models 
cover log(g)=4-4.5 and Te// = 900-lOOOK. The range 
in log(g) for HR 8799c and d are similar to that for b 
(log(g)=4.25-4.5 and 4-4.5), whereas their temperatures 
are slightly higher (1100-1200K and 1000~1200K). 

As illustrated by Figure [Tl] models with thick cloud 
layers provide far better fits to the SEDs of all three 
planets. Quantitatively, the x^ minima shrink by factors 
of 6, 2, and 5 for HR 8799b, c, and d compared to those 
for Model E fits. For HR 8799b and c, the minima ap- 
proach the formal 5-a confidence limit. For HR 8799d, 
multiple models have x^ minima less than the formal 3-a 
confidence limit. 

The righthand panels of Figure [Tl] illustrate why the 
thick cloud layer models are more accurate. For HR 
8799b, the best-fit models predict a flat, rising SED from 
1 to 1.5 /im, consistent with the planet's weak Y and J 
band emission. The best-flt models also predict stronger 
3.3 fim emission than in the standard model case and in 
better agreement with HR 8799b's measured [3.3] flux. 
While the best-fit model for HR 8799c underpredicts its 
J-band fiux while overpredicting its [3.3] and L' band 
flux, the discrepancies are weaker than in the standard 
cloud model case. With the exception of the CH4I filter 
data, which was not incorporated into our fitting, the 
best-flt thick cloud model (log(g) = 4.25, Tg// = HOOK) 
for HR 8799d accurately reproduces the planet's flux at 
every datapoint. 

Allowing the planet radii to vary by ± 10% slightly 
improves the model fits. More importantly, results in 
more models with x^ values below the formal 3(t and 5-tT 
confidence limits (Figure [T2|) . For these models, the HR 
8799b's range of best-fit models have log(g) = 4.25-4.5, 
and Te// = 900-lOOOK, and Cfc = 0.9-1.02; HR 8799c's 
have log(g)=4.25-4.5, Te// = 1100-1200K, and Ck = 
0.9-0.975; and HR 8799d's have log(g) = 3.75-4.5; Te// 
= 1000-1200K; and Cfc = 0.9-1.09. As before, the scal- 
ing factor for each model is correlated with the model's 
temperature compared to the cooling model estimates. 

3.3. Estimates for "Patchy " /Partly Cloudy Models 

The two models used to fit our data define limiting 
cases for the cloud structrure in planet atmospheres. 
The Model A thick cloud layer prescription fits the 
data for each planet far better. However, intermediate 
cases - with far thicker clouds than the Model E case 
but slightly thinner than Model A or a "patchy" cloud 
coverage - may be more physically re alistic. The two 
proces ses may be tied together: Acker nran and MarlevI 
()200H ) show that clouds may become patchy as they sed- 
iment below photospheric pressures. Near-IR photomet- 
ric variability detected from the T2.5 brown dwarf SIMP 
J013656. 5-1-093347 is consistent with grain free, cloudless 
regions and grain-bearing cloudy regions rotating in and 
out of view (jArtigau et al.l[2009l ). Cloud patchiness may 
also be important for d efining the L/T dwarf transition 
(e.g. iMarlev et aLll2010l and references therein). 
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We leave a detailed construction of such models to a 
future paper (Madhusudhan et al., in prep.) but here 
we qualitatively explore how intermediate cases may af- 
fect the predicted p lanet spectrum ( see also lMarlev et al.1 
12010). Similar to Burgasser et al. (2002), we follow a 
highly simplified, crude approach by combining weighted 
sums of Model A and E cloud prescriptions to approx- 
imate an atmosphere whose cloud thickness varies over 
the seeing disk of the planet^^. For simplicity, we com- 
pare two parameterizations: a "partly cloudy" approxi- 
mation where we weight the thick cloud model by 60% 
and a "mostly cloudy" approximation where we weight 
90% of the surface by the Model A case. 

Figure [T3] shows modeling results for these two cases 
compared against the thick cloud layer results for log(g) 
= 4 and 4.5. Our approximations yield smaller ^ min- 
ima for HR 8799b and c; models with partly/mostly 
cloudy approximations have the smallest . The best-fit 
model for HR 8799b has Tg// = 900K, consistent with 
the thick cloud layer model, while temperatures for HR 
8799c and d are lower by lOOK. 

While our approach is entirely ad hoc, it indicates 
that slightly weakening clouds compared to the limiting 
Model A case may provide better fits, at least for low 
surface gravity models (log(g) — 4). Madhusudhan et al. 
(2011) present a set of new atmosphere models with a 
range of cloud coverages intermediate between the Model 
A and E cases to explore how varying the cloud strength 
between these two extremes affects fits to the data. 

4. DYNAMICAL STABILITY ANALYSIS 

As shown bv iFabrvckv and Murray-Clavl ()2010l ) and 
iMoro-Martin et all ()2010D . stabilitv analvsis of the HR 
8799 system constrains the planet masses independently 
of planet cooling and atmospheric modeling. Here, we 
investigate the plausible mass range of companions im- 
posed by dynamical stability. Later, we will combine the 
results of these simulations with the implied mass range 
from atmospheric modeling to identify planet masses 
consistent with both atmospheric modeling and dynam- 
ical stability analysis. 

4.1. Procedure. 

Using the Swifter N-body code, an up dated version 
of the Swift package ([Duncan et al.l ll998\ we integrate 
the equations of motion for the HR 8799 planets. We 
adopt the Burlirsh-Stoer method to treat close encoun- 
ters. For all simulations we use an accuracy parameter 
of 10~^^ and dynamically evolve the system until one or 
more planets are ejected or un til 100 Myr is reached. 

We expand the analysis of iMaro is et al.l ()2011l) who 
searched for stable orbital configurations for two sets of 
planet masses incorporating HR 8799e - 5, 7, 7, and 7 
Mj; 7, 10, 10, 10 Mj for b, c, d, and e. We assume a 
single-2:l resonance between c and d for the former and a 
double-2:l resonance for d-c and c-b pairs for the latter. 
We hereafter refer to these sets of initial conditions as 

Technically, this is not physically realistic as the temperature- 
pressure profiles for cloud layer and cloud deck regions would 
bo discontinuous. On the other hand, for a giv en Tej j self- 
consi stent models with intermediate cloudiness ijMarlev et al.l 
[20T0I) have color -magnitude positions intermediate between the two 
extremes, broadly co nsistent with simple parameterizations (e.g. 
IBurgasser et al.ll200^ '). 



Cases A and B. To expand upon the lMarois et al.l (|2011[ ) 
investigation, we considered a wider range of masses for 
HR 8799bcde - 10, 13, 13, and 13 Mj - with the same 
double resonance configuration as Case B. We refer to 
this set of initial conditions as Case C. In all cases, we 
simply require the system to be stable for 30 Myr - the 
minimum age of HR 8799 - to be consistent with the 
data. 

We do two sets of 8000 simulations for each case. In the 
first set, we allow HR 8799e to vary in separation from 
13.1 AU to 15.7 AU. This allows us to identify general 
trends in the time to instability vs. separation for HR 
8799e. In the second set, we more finely sample initial 
orbital properties for the planets assuming a range of 14- 
15 AU for HR 8799e to better identify stable solutions. 

4.2. Results 

Figure [TJ] illustrates our simulation results. The top 
panel displays the time to instability for Case A. The 
bottom-left and bottom-right panels show the same plot 
for Cases B and C, respectively. The first set of simula- 
tions allowing HR 8799e to range from 13.1 AU to 15.7 
AU are shown as orange lines; the second set are shown 
as black lines. 

Our results show that the HR 8799 companions must 
have masses below the deuterium-burning limit based on 
dynamics alone. Case C configurations are typically only 
stable for 0.01 Myr and never stable for more than 10-20 
Myr. Because HR 8799 is a main sequence star, it cannot 
be as young as 10-20 Myr. Therefore, companion masses 
for HR 8799cde > 13 Mj and a mass for HR 8799b > 10 
Mj can be ruled out. 

Lower planet masses are strongly preferr ed on dynam- 
ical gr ounds, consistent with the results of iMarois et al.l 
(|20ll . Only seven Case B configurations are stable for 

30 Myr, nearly all of which require separations for HR 
8799e more than 1-cr different from the position implied 
by our astrometry. Only two are stable for 100 Myr, and 
these likewise require anomalously small separations. On 
the other hand, sixty Case A configurations are stable 
for 30 Myr. Three are stable for 100 Myr, one of which 
places HR 8799e at a separation consistent with our as- 
trometry. Our dynamical stability results are in rough 
agreement with iMarois et all (|2011t )'s results. They find 
12 solutions out of 10^ possible solutions stable for more 
than 100 Myr, where HR 8799e varies between 14.35 AU 
and 14.56 AU. We find 3 out of 1.6x10^ solutions are 
stable for 100 Myr over this semimajor axis range. 

In summary, we can rule out companion masses greater 
than 10 Mj for HR 8799b and 13 Mj for the others. 
The companions cannot be brown dwarfs. Systems with 
masses of 5 Mj for HR 8799b and 7 Mj for the others 
are characteristically far more stable than systems with 
larger masses. We fail to find any stable configuration 
with 7, 10, 10, and 10 Mj for HR 8799bcde's masses 
that place HR 8799e at a position consistent with our 
astometry. While our investigation is not exhaustive, it 
implies that masses of less than 7 Mj for HR 8799b and 
less than 10 Mj for the others are most plausible. 

5. DISCUSSION 

Our primary result in this paper is that the atmo- 
spheres of at least two and potentially all of the HR 8799 
planets do not easily fit within the empirical IR color 
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sequence for L/T type brown dwarfs of similar temper- 
atures, nor can they be well fit by standard atmosphere 
models used to infer the properties of brown dwarfs. 
Adopting realistic assumptions about planet radii, all at- 
mosphere model fits to data for HR 8799b and c are far 
poorer than any meaningful threshold identifying models 
consistent with the data. The models primarily fail by 
underpredicting the 3.3 /im flux and badly overpredict- 
ing flux at 1-1.3 /im. 

Our analysis suggests that having "thicker" clouds - 
ones with larger vertical extents ~ is key to reproducing 
the planets' SEDs. Compared to cloud structures as- 
sumed in standard L/T dwarf atmosphere models, these 
clouds are more optically thick at a given Teff, so they 
are visible (in the photosphere) at a lower Teff even 
though the cloud base is located far below at much higher 
pressures. Adopting a thick cloud layer prescription, we 
succeed in identifying models for each planet that quan- 
titatively are good-fitting models. Moreover, the tem- 
peratures of these models are consistent with simpler, 
presumably more accurate cooling model estimates. 

5.1. Comparisons with Previous Studies of HR 8799 

The most direct comparison to this work is the recent 
analysis of the HR 8799b K -band sp ectrum and 1.1-4.1 
/im photometry from Bo wler et al.l ( ^ 2010) whos e niod - 
eling formalism we largely follow. iBowler et ahl ()2010f ) 
also finds difficulties in using standard atmosphere mod- 
els to fit HR 8799b's SEP and interpret its properties 
(see also iMarois et al1l2008l ). Likewise, they find that 
temperatures inferred from standard atmosphere mod- 
els disagree with cooling model predictions and that the 
former require unphysically small radii. 

Our results indicate that including Y/z band data only 
exacerbates the already serious disagreement between 
standard cloud deck model predictions and the planet's 
SED. Our analysis confirms iBowler et al.l ()201Cl[ )'s infer- 
ence that HR 8799b's atmosphere is exceptionally dusty 
compared to field brown dwarfs. Our results extend this 
inference, indicating that HR 8799c and, plausibly, HR 
8 799d are also d usty c ompared to field brown dwarfs. 

iJanson et al.l ()201Cll) noted that while standard atmo- 
sphere models - the COND models in their case ~ can 
reproduce the mean brightness of HR 8799c's L'-band 
spectrum they incorrectly predict the spectral slope from 
3.9 fim to 4.5 /im. They cite greater atmospheric dust 
absorption and, especially, non-equilibrium carbon chem- 
istry as feat ures that may br ing the models into better 
agreement. iHinz et al.l (l2010l ) argue that incorporating 
non-equilibrium chemistry is necessary to reproduce the 
mid-IR photometry of HR 87 99bcd since the che mical 
equilibrium models they use (jSaumon pre- 
dict M-band fiuxes larger than the upper limits they re- 
port. 

Non-equilibrium carbon chemistry has little ef- 
fect on the near-IR p ortion of the SED (e.g. 
iHubenv and Burrow£l2007l) . Thus, our analysis indicates 
that thicker clouds - and, by implication, stronger atmo- 
spheric dust absorption - are far more important than 
non-equilibrium chemistry in reproducing the HR 8799 
planet 1-5 fim SEDs. Nevertheless, the HR 8799 planet 
atmospheres are plausibly not in local chemical equilib- 
rium. Since departures from chemical equilibrium alter 
the spectral structure at 4-5 /im, non-equilibrium chem- 



istry incorporated into thick or "patchy" cloud models 
may yield better fits to 1-5 /im photometry and mid- 
IR spectroscopy of the planets. Higher signal-to-noise L' 
band spectra and detections/more stringent upper lim- 
its at M will better identify evidence of non-equilibrium 
chemistry in the planets' atmospheres. 

5.2. Comparisons with Planet Evolution Models and 
Implied Masses 

Within the context of the iBurrows et al.l (|1997D planet 
cooling models, a particular combination of log(g) and 
Teff defines an object with a mass M and age t. Tak- 
ing the gravity and temperature range implied by our 
modeling at face value, we can then identify the mass 
and age range implied. Our modeling efforts succeed in 
yielding planets with physically realistic radii. However, 
if our range of log(g) and Teff were to imply wildly dis- 
crepant masses compared to cooling model estimates and 
dynamical stability requirements or widely varying ages 
our analysis would have solved one problem only to cre- 
ate comparably serious ones. 

Here, we combine all modeling results to identify the 
range of best-fit parameters and implied parameters - 
mass and age - from atmosphere models that we con- 
sider. We then determine whether the atmospheric and 
dynamical modeling constraints are consistent and, if so, 
what mass and age range they imply. 

• HR 8799b - The minimum value for HR 8799b 
for thick cloud models is 27.6 if w e allow the ra- 
dius t o vary by up to 10% from the IBurrows et al.1 
(|1997l ) values and 48.9 if we don't. For the 
"patchy" cloud approximation, the corresponding 

minima are 20.6 and 51.4. Considering the 
best-fit models passing the Ax^ threshold in each 
case, this range covers log(g) = 4-4.5 and Teff = 
800-lOOOK. Thus, our modeling yie lds log(g) = 4- 
4.5, T eff = 800-lOOOK. Using the IBurrows et al.l 
()1997l ) evolutionary models, this implies a mass and 
age range of M, t = 5 Mj, 30 Myr to 15 Mj, 300 
Myr. 

• HR 8799c, d, and e - The minimum values here 
for thick cloud models are 43.5 and 60.7 for c and 
5.7 and 5.3 with and without radius rescaling. For 
the "patchy" cloud approximation, the correspond- 
ing minima are 14-14.1 for c and 2.8-7.4 for 
d. For HR 8799c, the range of models passing the 
Ax^ threshold for the thick and patchy cloud pre- 
scriptions cover log(g) = 4-4.5 and Teff = lOOOK- 
1200K. This yields a mass/age range of 7 Mj, 30 
Myr to 15-17.5 Mj at 150-300 Myr. For HR 8799d, 
the range is log(g) = 3.75-4.5, Teff = 1000-1200K, 
yielding 5 Mj at 10 Myr to 15-17.5 Mj at 150-300 
Myr. Since HR 8799e likely has a bolometric lumi- 
nosity and K-L colors comparable to HR 8799c and 
d, its range of masses is plausibly consistent with 
those derived for HR 8799c and d. 

Dynamical constraints require that HR 8799b is less 
than 7 Mj and H R 8799cde are less than 10 Mj (Sec- 
tion 4 of this work; Ma rois et al1l2011[ ). The 5 Mj mass 
estimate for HR 8799d can be ruled out because the pri- 
mary star is on the main sequence and thus cannot be 
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10 Myr old. Coupled with the range in surface gravities 
and temperatures, the implied range in masses are then 
6-7 Mj for HR 8799b, 7-10 Mj for HR 8799c, and 7-10 
Mj for HR 8799 d. If HR 8799e's atmospheric proper- 
ties mirror those of c and d, its plausible range of masses 
is also 7-10 Mj. Conversely, for these ranges of masses, 
the surface gravities of HR 8799bcde should be no greater 
than log(g) « 4.25. 

These esti mates are cons i stent with c ooling model es- 
timates from iMarois et all ()2008l 1201 It ). For the lower 
end of the mass ranges, the system age corresponding to 
these models is « 30 Myr and puts HR 879 9's age on the 
low end of the 30-160 Myr range quoted bv lMarois et al.l 
(|2008). The (disfavored) high end of the mass range cor- 
responds to ^ 100 Myr-old objects. 

Despite our success in arriving at self-consistent an- 
swers for the planets' masses and ages, we strongly cau- 
tion against overinterpreting these results. Our results 
do not prove that, above the cloud base, the vertical 
density/pressure profile of clouds follows that of the gas 
as a whole (e.g. Si =0), as opposed to being truncated 
at higher pressures. Neither do our results prove that 
other models with slightly different assumptions about 
the clouds, grain particles, atmospheric chemistry, etc. 
provide better fits to the data. In particular, slight mod- 
ifications to our models may improve the fit at L' band, 
the datapoint responsible for much of the contribution 
for HR 8799b. Even within the context of our adopted 
physical models, our sampling in temperature and grav- 
ity is also too coarse to precisely estimate best-fit atmo- 
sphere parameters. 

On the other hand, our analysis provides compelling 
evidence for thick clouds, motivates future modeling 
work to test how different assumptions about thick 
clouds affect model fits to planetary atmospheres, and en- 
courages further observations of substellar objects to test 
these models. Madhusudhan et al. (2011) will develop 
and better assess model fits for varying cloud strengths 
and more precisely and accurately determine tempera- 
tures and gravities for the HR 8799 planets and other 
planetary-mass objects. 

5.3. Constraints On The Formation of the HR 8799 
Planetary System 

The planets' large masses and wide orbits make them 
a particularly interesting probe of planet formation. 
The favored theory invoked to explain the fo rmation 
of ga s giant plane t s is c ore accretion (e.g. iMizunol 
19801 iPollack et al.l Il996t iKenvon and BromlevI 120091 
Chambers et al.l|2010[ ). where cores that have grown to 
« 5-10 rapidly accrete much more massive gaseous 
envelopes . Alternati vely, planets could form by disk in- 
stability (jBossI [T997I and later papers), where the pro- 
toplanetary disk is massive and gravitationally unstable, 
forming multiple self-gravitating clumps of gas that coa- 
lesce into bound, planet-mass objects. 

HR 8799's planets are often des cribed as confound ing 
either planet formation model fe.g. lMarois et al]l2011t ) or 
being clear e xamples of disk instability-form ed planets, 
as claimed bv iDodson- Robinson et al.l (|2009D . They find 
that cores at distances characterizing the HR 8799 plan- 
ets cannot reach 10 M^ in mass to undergo runaway 
gas accretion even under the most favorable conditions. 
They claim that planet-planet scattering cannot create 



stable, wide-orbit systems like HR 8799's. They conclude 
that massive, wide-separation gas giants like HR 8799bcd 
form by disk instability and "can certainly rule out core 
accretion" . 

Critical to iDodson - Robins on et all 's conclusion is their 
treatment of the core growth rate. The growth rate 
strongly depends upon the planetsimal approach veloc- 
ity, which they fix at v^ = niihui- They claim this ve- 
locity yields an "optimistically high" growth rate. Their 
formalism implicitly assumes that planetesimals have an 
isotropic velocity dispersion (v^ ~ v^), which is valid 
as long as the scale height of planetesimals accreted by 
cores (vz/rj) is larger than the core's impact parameter, 
Rcore^/il + 0) (|Rafikovl[200l . where 9 is the Safranov 
number. However, if the planetesimals are dynamically 
cold such that v^ < ^nRH ui (where p = RcorefRHiu) , 
this conditio n is violated (jPones and Tremaind 119931 : 
lRafikovl[200l . The core can then accrete the entire ver- 
tical column of planetesimals at a vastly higher r ate since 
accre tion is now essentially two-dimensional (|Rafikovl 
[2OOI . 

As a result. [Podson- Robinson et al.l (|2009[ ) catastroph- 
ically underestimate the maximum growth rate by a 
factor of p~^/^, or up to 114, 85, and 68 at the po- 
sitions of HR 8799b, c, and d (cf. Equations 78, 
80, and 82 in Rafikov 2004; see also Rafikov 2010)^6. 
Detailed numerical simulations confirm that this rapid 
growth phase can be reached if col lisional fra gmentation 
and g as dra g are properly treated (iKenvon and Bromle^ 
120091) . The Dodson-Robinson et al.l planet-planet scat- 
tering simulations also were conducted assuming gas free, 
planetesimal-free conditions and assumed that planets 
could not further grow after scattering. However, gas 
drag and dynamical friction from planetesimals are criti- 
cally important as they promote orbit c ircularization and 
stability (e.g. iGoldreich et al.l 12004 iFord and Chiang! 
^OOTh ". Cores with masses sufficient for rapid gas ac- 
cretion can circularize after being scattered to the outer 
disk (Bromley and Kcnyon 2011, S. Kenyon 2010, pvt. 
comm.). Simulations by Thommes et al. (in prep.) show 
that the HR 8799 planet cores could acquire most of their 
gas after scattering. 

The mass ratio and semimajor axis distribution of 
wide planets and low-mass brown dwarfs may help con- 
stra in the format i on m echanism for HR 8799's plan- 
ets (jKratter et al.l[2bl0t ). Core accretion preferentially 
forms planets with smaller masses and orbital sepa- 
rations, while disk instability has difficulty producing 
lower-mas s gas g i ants a nd forming them close to the 
star (e.g. iRafikovl 120051 : iKratter et'ali I2OIOI ). There- 
fore, if HR 8799bcde formed by core accretion (disk 

IS At first glance, Equation (17) in IRafiko vl 1 I2OIOI) appears 
to imply that the limiting distance for cor e accreti on in s hear- 
dominated growth is comparable to [Dodson-Robinson et al,Ps es- 
timate (44 AU vs. their 20—35 AU). However, Rafikov's result of 
44 AU is valid for a Minimum Mass Solar Nebula case ( Havashil 
|T981). Adopting i nitial a ssumptions more comparable to those that 
IDod son-Robinson et al] assumes - e.g. a disk more massive than 
the Minimum Mass Solar Nebula or a longer-lived one with T^j^afc = 
5 Myr instead of 3 Myr- implies that gas giants can in some cases 
form by core accretion at separations comparable to HR 8799c and 
b. 

In fairness, they clearly acknowledge that their study does 
not consider planet-planet scattering in a gaseous disk, which may 
result in a more favorable outcome for core accretion. 
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instability), they should comprise the high-mass ex- 
trema (low-mass tail) of a population continuous with 
radial-velocity detected planets (brown dwarf compan- 
ions). Using our new re sults for the m asses of the HR 
8799 planets, we update 'Kratte r et all s plot comparing 
planet and brown dwarf distributions. We also add the 
planet-mass c ompanions to 1RX JS1609. 1-210524, and 
2M J044144b (jTodorov et al.llMol 5-10 Mj, 15 AU); the 
planet/brown dwarf companion to GSC 06214-00210B 
(14 Mj, ~ 300 AU 'Irelan d et all |2010D : and the low- 
ma ss brown dwarf com panion GJ 758B (25-40 Mj, 44 
AU lCurrie eTal\\2()T(h . 

As shown by Figure [151 the revised masses for the HR 
8799 planets and the addition of HR 8799e expand the 
space between them and brown dwarf companions (aster- 
isks). Visually, they join with the distribution of closer- 
separation planets plausibly formed by core accretion. 
The other new companions (red triangles) are continuous 
with brown dwarfs that may form by disk fragmentation. 

While core accretion - especially when coupled with 
planet-planet scattering - may form the HR 8799 plan- 
etary system, HR 8799-like systems are still plausibly 
uncommon. The Gemini Deep Planetary Survey of 85 
nearby, young (mostly solar-mass) stars was typically 
sensitive to 2 Mj plan ets at 40-200 AU yet failed to 
detect any pLafrcnicr e et al.l [200 7b'). Similarly, non- 
detections fr o m th e deep (M < 1 Mj) survey from 
iKasper et al.l (|2007[ ) showed that the giant planet pop- 
ulations detected at small separations (a < 4 AU) by 
RV surveys cannot extend to separations larger than ^ 
30 AU. More massive stars like HR 8799 likely have 
more massive disks, which aid gas giant planet forma- 
tion. However, their disk s also dissipate much more 
rapidly ([Currie et al.ll2009l) : even if critical core masses 
are reached, the leftover mass of gas may be small. More- 
over, rapid core growth results from being able to frag- 
ment and then dynamical cool the surrounding planetesi- 
mal population. The current state-of-the-art simulations 
sho w that this requires Pluto-m ass cores to start with 
(e.g. lKenvon and Bromlevll2009[ ). yet the formation time 
for Pluto-mass ob jects at wide separations may be long 
(e.g. lRafikovll2010[ ). Thus, forming HR 8799-like systems 
by core accretion is difficult, though not impossible, and 
probably happens infrequently. 

5.4. Implications for the Atmospheres of Other 
Substellar Companions: A Possible Fundamental 
Difference Between Planetary- Mass Objects and M 
> 15-20 Mj Brown Dwarfs 

In some ways, the difficulty in reproducing the IR 
SEDs of the HR 8799 planets mirror difficulties in mod- 
eling other planetary-mass objects detected prior to HR 
8799bcde. In particular, 2M 1207b also appears dis- 
crepant c ompared to sta i idard atmosphere models as 
noted in iMohantv et al.l ([20071 ) and discussed in this 
work. Like HR 8799b, 2M 1207b is no ticeably under- 
lumin ous (~ 2.5 mags) in the near-IR ([Mohantv et al.l 
120071 this work). 

To explain 2M 1207b's properties, IMohantv et al.l 
(|200l propose that the object is occulted by an edge- 
on disk with large, gray d ust grains. Alternatively, 
iMamaiek and Meved ([2007) propose that 2M 1207b's 
properties can be explained as resulting from a re- 
cent protoplanet-protoplanet collision. Comparing high- 



resolution spectra o f 2M 1207b t o the DUSTY atmo- 
sphere models from lAllard et al.l ([20011 ) . [Patience et al.l 
([2010() identify a pro blem similar to that noted for mod- 
eling HR 8799b from [Bowler et al.l ([2O1O0 and this work. 
Namely, allowing the object radius to freely vary yields 
best-fit radii far smaller (fa 0.5 Rj) than is physically 
plausible (cf. [Burrows et al.lll997D . [Patience et all ([2010[ ) 
also conclude that extinction from an edge-on disk com- 
prised of gray dust grains is also a viable scenario. 

For the same reasons - underluminosity/red colors - 
a disk origin also has been proposed to explain the IR 
SED of HR 8799b and (to a lesser exten t ) c a nd d 
([Marois et al.ll2008l ). However. [Marois et al.l ([2008[ ) con- 
sider the chance alignment of an edge-on circumplanc- 
tary disk to be unlikely, especially given that the sys- 
tem is viewed nearly face on. Even more unlikely is the 
chance that circumplanetary disks or recent protoplanet 
collisions explain the near-IR properties of two to four 
separate planets in two systems with very different ages 
and primary star properties. 

Given the success in better reproducing HR 8799bcd's 
SEDs with thick cloud models and the similarity between 
HR 8799b and 2M 1207b, it is more plausible that the 
latter's near-IR spectrum is likewise explained by thick 
clouds. If this is generally true of planetary-mass ob- 
jects, thicker clouds may constitute the primary differ- 
ence between the atmospheres of massive planets and 
brown dwarfs, at least over the gravity and temperature 
range enclosed by the HR 8799 planets and 2M 1207b 
(e.g. log(g) = 3.75-4.5, Te// =900-1600K). Since thicker 
clouds affect the color-magnitude positions of substellar 
objects it is quite possible the Model A 'thick cloud' se- 
quence extending to HR 8799b and 2M 1207b from the 
nominal L/T dwarf boundary continues on to even cooler 
temperatures (e.g. Te// ~ 700-900K). Since thick clouds 
present reshape the spectral structure at ~ 1.6 fim (e.g. 
in the methane band), they may also affect the L/T 
dwarf transition, which is alr eady known to be depen- 
dent upon surface gravity (e.g. iMetchev and Hillenbrandl 
[20MlLuhman et al.l[2007l ). 

5.5. Future Work 

Our study motivates the development of a suite of new 
atmosphere models with clouds intermediate in thick- 
ness between the Model E cloud deck and Model A 
thick cloud layer prescriptions. Adopting these mod- 
els as fiducial models, we can revisit the (secondary) ef- 
fects of surface gravity, metallicity, and non-equilibrium 
chemistry on the atmospheres of planetary mass objects, 
complementing similar investigations for brown dwarfs 
(Allard et al. 2001; Marlcv ct al. 2002; Burrows ct a]] 
l2006t [Hubenv and Burrowsl |2007|) . These models will 
be developed and applied to HR 8799bcde and other 
planetary-mass objects in upcoming papers (Madhusud- 
han et al. 2011, in prep.) and may provide a use- 
ful comparison to planet parameters derived from cool- 
ing models (e g. 'Burro ws et "allll997[: [Baraffe et al]|2003[: 
iFortnev et aI[[2007L ,2001 1 

New observations at 1-5 fim will provide better con- 
straints on the HR 8799 planet atmospheres. In addi- 
tion to more sensitive data at Y band and [3.3], Fig- 
ure [H (lower-left panel) implies that thick-cloud atmo- 
spheres may have far stronger emission at ^ 2.3 /im 
and 3.0 fim than standard models predict. This wave- 
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length range can be probed for at least HR 8799bcd 
by current ground-based facilities such as VLT/NaCo, 
Keck/NIRC2, and MMT/Clio. Integ ral field spectro- 
graph s on Gemini Pl anet Imager (GPI Macintosh et al] 
120081 ) and SPHERE (Beuzit ct al. 2008') will sample the 
1-2.5 /im SED region with exceptional sensitivity and 
thus provide a detailed comparison between observed and 
predicted atmospheric properties of all planets. 

Finally, ongoing collaborations such as the IDPS 
survey (Marois et al., in progress) and Gemini/NICI 
(|Liu et al.ll2010[ ) will better probe the frequency of wide, 
massive (~ 5-13 Mj, > 30 AU) around nearby stars. GPI 
and SPHERE will probe 1-5 Mj planets at even smaller 
separations (e.g. 5-30 AU). These surveys will produce 
a far more complete census of Jupiter-mass planets to 
better determine their ubiquity and constrain how the 
formation of planets like HR 8799's compare to that ex- 
pected for lower-mass planets at smaller separations and 
wide-separation, low-mass brown dwarfs. 
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TABLE 1 
Observations 



Telescope/Instrument 


Filter 


Date 


Exposure Time (s) 


Field Rotation (degrees) 


Detections 


MMT/CLIO 


L' 


November 21, 2008 


5460 


83.88 


b, c, d 




M 


November 21, 2008 


9600 


31.8 






[3.3] 


September 12, 2009 


6780 


128.53 


bi, c 


Subaru/IRCS 


z 


August 15, 2009 


4200 


172 


bi 




J 


August 15, 2009 


1080 


7.4 


b 


VLT/NaCo 


K, 


October 8, 2009 


6185 


63 


b, c, d, e 



Note. — Note (1) - Companion has a low signal-to-noise detection because it is intrinsically faint 
(photon- noise dominated region). 



TABLE 2 

ASTROMETRY (E["], N["] POSITION) 



Date 11-21-2008 8-15-2009 9-12-2009 10-8-2009 
Planet 

HR 8799b 1.532 ± 0.02, 0.796 ± 0.02 1.536 ± 0.01. 0.785 ± 0.01 1.538 ± 0.03, 0.777 ± 0.03 1.532 ± 0.007, 0.783 ± 0.007 

HR 8799c -0.654 ± 0.02, 0.700 ± 0.02 - -0.634 ± 0.03, 0.697 ± 0.03 -0.627 ± 0.007, 0.716 ± 0.007 

HR 8799d -0.217 ± 0.02, -0.608 ± 0.02 - - -0.241 ± 0.007, -0.586 ± 0.007 

HR 8799e - - - -0.306 ± 0.007, -0.217 ± 0.007 



Note. — The 8-15-2009 astromctry listed for HR 8799b comes from the .J band data because this data yields a higher signal-to- noise detection. 



TABLE 3 
Photometry 



Filter 


z 


J 


H 


CH45 


CH4i, 




[3.3] 


L' 


M 


A {fJ.m) 


1.03 


1.248 


1.633 


1.592 


1.681 


2.146 


3.3 


3.776 


4.8 


Planet 




















b 


18.24 ± 0.29 


16.52 ± 0.14 


14.87 ± 0.17 


15.18 ± 0.17 


14.89 ± 0.18 


14.05 ± 0.08 


13.96 ± 0.28 


12.68 ± 0.12 


> 11.37 


c 


> 16.48 


14.65 ± 0.17 


13.93 ± 0.17 


14.25 ± 0.19 


13.90 ± 0.19 


13.13 ± 0.08 


12.64 ± 0.20 


11.83 ± 0.07 


> 11.22 


d 


> 15.03 


15.26 ± 0.43 


13.86 ± 0.22 


14.03 ± 0.30 


14.57 ± 0.23 


13.11 ± 0.12 


> 11.63 


11.50 ± 0.12 


> 11.15 


e 












12.89 ± 0.26 




11.61 ± 0.12 





Note. — Magnitudes list ed are the absolute magnitude of the companions, assuming a distance of 39. 4 pc. (1) H, CH4S, CH4L, and K band photometry 
for HR 8799bcd taken fr om JMarois et alj200gH . J band photometry for HR 8799c and d also taken from lMarois et al.1 l |200gr . L' band (3.8 fim) photometry 
for HR S799e taken from fMaroi^^t^!n '' j20 111) . Photometrv/upper limits at 3.3 fim, L' band and M band (4.8 fim) for HR 8799bed are taken from this work. 



TABLE 4 

Adopted Photometry for Other Planet-Mass Objects and Low-Mass Brown Dwarfs 



Companion 


D (pc) 


J 






H 


K 






L' 


References 


2M 1207b 


52.4 


16.40 


± 


0.2 


14.49 ± 0.21 


13.33 


± 


0.11 


11.68 ± 0.14 


1,2 


1RXJ1609.1-210524 


140 


12.17 


± 


0.12 


11.139 ± 0.07 


10.44 


± 


0.18 


9.14 ± 0.3 


3,4 


AB Pic b 


47.3 


12.80 


± 


0.10 


11.31 ± 0.10 


10.76 


± 


0.08 


-99 


5 


HD 203030b 


40.8 


15.08 


± 


0.55 


13.80 ± 0.12 


13.16 


± 


0.10 


-99 


6 



Note. — All magnitudes listed arc a bsolute magni tude s. References - 1) [ Chauvin et al.V2004"), 
2) iMohantv et al . (2007), 3) Lafrcnic re et al.l (|2008l '). 4) [Ireland et al.l (|20ia ), 5) .Chauvin ct al 
(|2005t ) ■ 6) iMetchev and HiUenbrandl (|2006l ) . 



TABLE 5 

Standard Model (Model E) Fitting Results 



Model Run 


¥^ ■ 


log(g), T,ff (for xLn) 






log(g), log(Te//) (x^ < Ax^) 


HR 8799b 












Model E soleir, 3x 


279.0 


4.5, 900K 


1 


291.85 


4.5-5, 900-lOOOK 


Model E solar, 3x 


264.5 


4.5, 900K 


0.91 


286.4 


4.5, 900K; 5, lOOOK 


(Cfc = 0.9-1.1) 












Model E solar, 3x 


36.6 


4.5, 1400 


0.34 


58.5 


4, 1400K; 4.5, 1300-1500K; 5, 1400K 


(Cfc =0.2-2) 










HR 8799c 












Model E solar, 3x 


120.8 


5, 1200K 


1 


142.7 


4.5, HOOK; 5, 1200K-1300K 




71.0 


O.U, -LOUUxv 


0.9 


92.5 


O, -LOUUlv 


(Cfc=0.9-l.l) 












Model E solar, 3x 


17.6 


4.5, 1400 


0.54 


39.5 


4, 1300-1700K; 4.5, 1300-1700K; 5, 1400-1700K 


(Cfc=0.2-2) 








HR 8799d 












Model E solar, 3x 


17.0 


4.5, HOOK 


1 


38.5 


4.5, HOOK; 5, 1200K 


Model E solar, 3x 


17.0 


4.5, HOOK 


0.9 


38.5 


4.5, 1100-1200K; 5, 1200-1300K 


(Cfc=0.9-l.l) 












Model E solar, 3x 


10.91 


4.5, 1300K 


0.64 


32.8 


4, 1300-1600K; 4.5, 900— 1700K; 5, 1200-1700K 


(Cfc=0.2-2) 









Note. — Where the metal rich models are considered (first three rows for each planet), they always provide the smallest 
value. 



TABLE 6 

Thick Cloud Model (Model A) and "Patchy" Cloud Approximation Fitting Results 
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Model Run 



log(g), T,ff (for xLn) Cfe log(g), log(T,/^) (x^ < Ax^) 



HR 8799b 

Model A solcir 
Model A solar 
(Cfe=0.9-l.l) 

Model A/E solar 

(60% thick clouds) 

Model A/E solar 

(60% thick clouds, Cfc=0.9-l.l) 

Model A/E solar 

(90% thick clouds) 

Model A/E solar 

(90% thick clouds, 0^=0.9-1.1) 

HR 8799c 

Model A solar 

Model A solar 
(Cfc=0.9-l.l) 

Model A solar 
(60% thick clouds) 
Model A solar 

(60% thick clouds, Cfc=0.9-l.l) 

Model A solar 
(90% thick clouds) 
Model A solar 

(90% thick clouds, Cfc=0.9-l.l) 

HR 8799d 

Model A solar 

Model A solar 
(Cfc=0.9-l.l) 

Model A solar 
(60% thick clouds) 
Model A solar 

(60% thick clouds, Cfe=0.9-l.l) 
Model A solar 
(90% thick clouds) 
Model A solar 

(90% thick clouds, Cfc=0.9-l.l) 



48.9 4.25, 900K 
27.6 4.5, lOOOK 



91.5 4, 900K 

85.5 4, 900K 

51.4 4, 900K 

20.6 4, 900K 

60.7 4.5, HOOK 

43.5 4.5 1200K 

14.1 4, lOOOK 

14.0 4, lOOOK 

51.6 4, lOOOK 
48.0 4, HOOK 



5.7 
5.3 

7.4 
2.8 
8.8 
3.1 



4.25, HOOK 
4.5, 1200K 

4, lOOOK 
4, lOOOK 
4, lOOOK 
4, lOOOK 



1 70.8 4, 900K; 4.25, 900K; 4.5, 900-1000K 
0.9 47.7 4, 900K; 4.25, 900K; 4.5, 900-lOOOK 



1 111.6 4, 800K 

1 105.6 4, 800-900K 

1 71.5 4, 800-900K 

1 40.7 4, 800-900K 



1 82.6 4.25, HOOK; 4.25, HOOK 

4.5, 1100K-1200K 
0.9 63.6 4.25, HOOK; 4.5, 1100-1200K 



1 34.5 4, lOOOK 

1 34.1 4, 1000- HOOK 

1 71.7 4, lOOOK 

1 68.1 4, lOOO-llOOK 



1 25.8 3.75, lOOOK; 4, lOOO-llOOK 
4.25, HOOK; 4.5, 1100-1200K 

0.98 27.2 3.75-4, lOOO-llOOK; 4.25, 1000-1200K 
4.5, 1100-1200K 

1 27.5 4, lOOO-llOOK 

1 22.9 4, lOOO-llOOK 

1 28.9 4, lOOO-llOOK 

1 23.2 4, lOOO-llOOK 
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Fig. 1. — VLT/NaCo image of the HR 8799 planetary system. Previously detected planets — HR 8799b, c, and d — are easily visible at 
high signal-to-noise. At ~ 0.375" separation, we detect an add itional object consistent with being a fourth planet orbiting HR 8799 - HR 
8799e. This same object was independently detected bv lMarois et al.l II2011I ) and confirmed to be a fourth planet. HR 8799e (and d, to a 
lesser extent) appear slightly smaller than b and c because of point source self subtraction inherent in LOCI processing. 
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Fig. 2. — Astrometry for HR 8799e compari ng positions from [Marois et al.1 pOlU ') and from our work. The arrow identifies the direction 
to the HR 8799 primary. The two points from lMarois et ahl l l201ll ) nearest to our October 8, 2009 me asurement were taken in August 2009 
and November 2009, respectively. Our astrometry are consistent with those from lMarois et aP 11201 ll ) within errors. 
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Fig. 3. — Subaru/IRCS images (left panels) and signal-to-noise maps (right panels) obtained at J band (top panels) and z band (bottom 
panels). The poor field rotation and short integration time in J limit our detection to HR 8799b. Despite over an hour of integration time 
in z band, we marginally detect HR 8799b but fail to detect the other planets. 




Fig. 4. — MMT/Clio image of the HR 8799 system at L'/3.8 iim. The three planets are clearly visible and all are detected at SNR > 5. 



18 





Fig. 5. — MMT/Clio image at 3.3 fim. (Left panel) The image shown with a high dynamic range to more clearly show the detection of 
HR 8799c. (Right panel) The image with a smaller dynamic range to better illustrate the marginal detection of HR 8799b. 





Fig. 6. — MMT/Clio image at M'/4.8 fim shown with slightly different procedures for PSF subtraction (left panel, simple ADI reduction; 
right panel, radial profile subtraction with LOCI reduction) and different dynamic ranges (left panel, high dynamic range; right panel, low 
dynamic range to show residual noise). The circle identifies the centroid position of HIl 8799c in the L' image obtained on the same night. 
While a local peak appears near the position of HR 8799c, we do not identify any > 3-cr peaks consistent with any of the planets in these 
images. 
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Fig. 7. — Plot of the point source self subtraction for LOCI as a function of separation for our MMT/Clio data for a simple ADI reduction 
and our LOCI reduction. LOCI attenuates more flux, especially interior to 0.75". 
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Fig. 8. — Color-magnitude diagrams comparing the HR 8799 planets with field L dwarfs (black dots) and T dwarfs (grey dots) and other 
planetary or very low-mass brown dwarf companions (squares). In Ks/Kj-L' (top-left panel), the planets follow the L/T dwarf sequence. In 
at least one of the diagrams including Y, J, and H-band data (top-right panel; bottom panels), the planets are red/underluminous compared 
to the empirical L/T dwarf sequence and the synthetic L/T dwarf colors from lBurrows et al, (200^) for a range of metallicities and gravities. 
The positions for other planetary- mass/low-mass brown dwarf companions also depart from the L/T dwarf sequence, especially 2M 1207b. 
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Fig. 9.— 



Comparing SEDs for different cloud prescriptions (no clouds, Model E, and Model A) at a given temperature and gravity. 
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Fig. 10. — Fitting results for the standard cloud deck models assuming the IBurrows et al.l (119971 ) radii. The top-left panel show the 
distribution of vs. T^f j for model fits to HR 8799b with a range of surface gravity and metallicity. The top-right panel compares the 
HR 8799b planet SED to the model with the smallest value. The bottom panels compare the HR 8799c (left) and HR 8799d (right) 
SEDs to the best-fit models for these data. In the SED comparisons, the horizontal magenta lines identify the flux of the model in the 
photometric filters convolved over the filter function. The width of the magenta line corresponds to the width of the filter. 
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Fig. 11. — Fitting results for the Model A, thick cloud layer prescription with a modal particle size of 60 fj,m. The lefthand panels show 
the distributions for each planet while the righthand panels compare the planet SEDs to the best-fit models in each case. Compared to 
the Model E, standard cloud deck fits, these models yielded smaller minima and better fits to the data. 
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Fig. 12. — Same as Figure [TT1 except allowing the planet radius to depart by ± 10% from the lBurrows et aP l|1997l ) values. 
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Fig. 13. — Fitting results for our simple approximation of a "patchy" cloud atmosphere. In all right-hand panels, the displayed best-fit 
SEDs have = 0.9. 
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Fig. 14. — The time to dynamical instability vs. semimajor axis of HR 8799e for three separate sets of masses and orbital resonances. 
In all plots, the orange lines denote our set of simulations allowing HR 8799e to vary between 13.1 AU and 15.7 AU, while the black lines 
denote our simulations that restrict HR 8799e to be between 14 and 15 AU while more finely sampling the range of other orbital parameters 
(e.g. mean anomaly, longitude of periastron). The horizontal dot with error bars identifies the l-cr range of projected separations for 
HR 8799e from our work. In Case A (top panel), HR 8799bcde have masses of 5, 7, 7, and 7 Mj. Case B (bottom-left) corresponds to 
planet masses of 7, 10, 10, and 10 Mj and Case C (bottom-right) correspond to 10, 13, 13, and 13 Mj. The density of the bars appears 
anomalously low for Case C because many simulations have instability times less than 10^ years. 
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Fig. 15. — Updated version of the mass ratio vs. orbital separation plot from iKratter et al.' f2010') incorporating our revised masses for 
HR 8799bcd (magenta dots). We also include HR 8799e, assigned a mass o f 7 M,/ fromiM arois ct al. (2011), and displayed as the left- 
most magenta dot. The [3 Pic planet is specifically identified as a green dot l ILagranee et al.ll201Cll ). Fomalhaut b is a downward-pointing 
magenta triangle (Kalas et al. 2008; Chiang ct al. 2009). Substellar companions discovered after or not included in the lKratter et aU 1 120101 ) 
pubUcation - 1RXJS1609. 1-210524, GJ 758B, 2M J044144b, and GSC 06214-00210B - are plotted as red triangles. Black crosses, purple 
crosses and purple squares denote radial velocity, transit, and microlensing-detected planets around stars with three mass bins: M* > 0.4 
Mq, Mq = 0.1-0.4 M0, and M* < 0. 1 M.p). F or direct comparisons and simplicity, wc plot same the population of exoplanets not detected 
by direct imaging as that used Kratt eFet al.l ([2010) (e.g. we do not include planets discove red by RV or tran sits since the publication 
of this paper). Black asterisks denote the sample of substellar companions listed in iZuckerman and Sond (120091) . Jupiter and Saturn are 
plotted as blue diamonds. 



